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The efTect of noise on the process of high-speed remagnetization of vortex state of a pentagonal 
array of five circular magnetic nanoparticles is studied by means of computer simulation of Landau- 
Lifshits model. The mean switching time and its standard deviation of the reversal between the 
counterclockwise and clockwise vorticities have been computed. It has been demonstrated that with 
the reversal by the pulse with sinusoidal shape, the optimal pulse duration exists, which minimizes 
both the mean switching time (MST) and the standard deviation (SD). Besides, both MST and 
SD significantly depend on the angle between the reversal magnetic field and pentagon edge, and 
the optimal angle roughly equals 10 degrees. Also, it is demonstrated that the optimization of the 
angle, duration and the amplitude of the driving field leads to significant decrease of both MST and 
SD. In particular, for the considered parameters, the MST can be decreased from 60 ns to 2-3 ns. 
Such a chain of magnetic nanoparticles can effectively be used as an element of magnetoresistive 
memory, and at the temperature 300K the stable operation of the element is observed up to rather 
small size of nanoparticles with the raduis of 20 nm. 



Studies on arrays of magnetic nanoparticles have 
arisen growing interest due to their potential applica- 
tions in spintronics devices and especially in the past 
decade with the development of modern technologies, 
allowing to produce regular arrays of magnetic nanopar- 
ticles with minimal sizes of order 20 nm. The peculiari- 
ties of the magnetic collective behavior have essentially 
been investigated by numerical calculations taking into 
account the relevant interparticle interactions. One of 
the directions of such investigations is the study of the 
magnetic flux closure (MFC) state or the magnetization 
vortex curling state realized within a single magnetic na- 
noelement. Such artificially designed arrays of magnetic 
nanoparticles can have clear application for magnetic 
logic devices^. On the other hand, one of perspec- 
tive realizations of the magnetic random access memory 
(MRAM) uses magnetic rings with vortex-like distribu- 
tion of magnetization. The interest to such vortex states 
is due to two reasons: first, the crosstalk fields are small 
in this state, allowing high density of elements. The 
second, the remagnetization of such elements can simul- 
taneously be driven by magnetic field of a current, hav- 
ing vortex-like structure, and also with the use of the 
spin torque effect. However, considering large circular 
nanoparticles^, it should be noted that it is practically 
impossible to perform their optimization in frequency re- 
sponse and stability to thermal noise due to large calcu- 
lation time. That is why it is important to study fluctu- 
ational effects in systems of only few magnetic nanopar- 
ticles where MFC states are realized. In Refii the cluster 
of five circular magnetic nanoparticles forming a ring (let 
us call it the pentagonal array) has been considered and 
the approach for the formation and manipulation of the 
chirality of a magnetization flux closure (MFC) state has 
been proposed. However, the study has been focused on 
the hysteretic loops and such interesting topic as switch- 
ing times between opposite vortex states has not been 
investigated. Besides, the effect of thermal fluctuations 
on the stability of MFC state and switching times were 
not discussed. 

Due to small size the nanoparticles can often be con- 
sidered as single-domain particle o^ since the inho- 



mogeneous distribution of the magnetization is energet- 
ically disadvantageous. Due to the complexity of the 
model, described by the time-dependent Landau-Lifshits 
equation with noisc^^, starting from the seminal paper- 
mostly the relaxation times of magnetization in the case 
of static magnetic fleld have been studiecU^^. Ref-si^ii^ 
were focused on the effect of noise and high-frequency 
rotating magnetic field on the magnetic dipole dynam- 
ics, while Refi^ was devoted to investigation of stochas- 
tic resonance effect. However, the practically interesting 
effect of fluctuations on the high-speed switching of mag- 
netic single-domain particles has been investigated only 
recently^ (2^. 

In the present paper the effect of thermal fluctuations 
on the process of high-speed remagnetization of vortex 
state of a chain of five circular uniformly magnetized 
nanomagnets, located in vertices of a pentagon, is stud- 
ied by means of computer simulation of Landau-Lifshits 
model. It is focused on the investigation of statistical 
characteristics of the reversal process with the aim to 
find an optimal regime of reversal between different vor- 
tex states with the smallest mean reversal time and the 
standard deviation. 

The dynamics of a magnetic nanoparticle of an array 
is described by the Landau-Lifshits equationi^: 
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where is the magnetization of a particle, l} is the 
effective magnetic field, 7 is the gyromagnetic constant, 



is the satu- 



a is the damping, /3 = 1 + 

ration magnetization. The effective magnetic field con- 
tains the following components: it — Hd + + Hr, 
where Hd is the field of magnetic-dipole interaction, 
He is the external field, and Ht - fluctuational field. 
The fluctuational field is assumed to be white Gaus- 
sian noise with zero mean and the correlation function: 
{H{t)nH{t')Tj) = ^jT^Sit - t')S,,, where k - Bohz- 
mann constant, T is the temperature, and V — |7r(Ai?)^ 
is the volume of the magnetic particle. Here, as ir*^, 
each magnetic nanoparticle is assumed to have a core- 



shell structure, with a soft ferromagnetic core encapsu- 
lated by a thin shell of nonmagnetic material. The core- 
shell structure is specified by two geometric parameters, 
which are the outer radius R and the core-to-particle 
radius ratio, A = Rcore/R, where Rcore is the radius of 
the magnetic core. Adjacent nanoparticles in the array 
are physically in contact to form a compact array struc- 
ture so that the outer diameter of each particle 2R is 
the particle spacing. Since the magnetic core is made of 
soft ferromagnetic material in perfect spherical shape, 
both the shape and the magnetocrystalline anisotropies 
are small and are not accounted for in the calculation. 
Owing to the presence of the nonmagnetic shell, the in- 
terparticle interaction is limited only to the long-range 
dipolar magnetic interaction. Below we will study two 
cases A = 0.9 and A — 1.0. The magnetic-dipole inter- 
action field for the i-th dipole has the following form: 
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where Vr = ^tt{\R)^ /R^, l^i is the radius- vector 
of i-th point, rij = I?*!— fj| is the distance be- 
tween positions i and j in units of the particle ra- 
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sin(^/jv) (sin ^ - sin ^) for A^=5 in our case. If the 
switching magnetic field is absent, due to magnetic- 
dipole interaction, the vorticity will be formed almost 
instantly, see Fig. [TJ Its direction depends on the sym- 
metry of the system of dipoles and it is possible to con- 
trol the sign of the vorticity using uniform external mag- 
netic field, exceeding the interaction field. For example, 
if the initial state is counterclockwise (see the left in- 
set of Fig. 1), and an external field is applied along 
y-axis and along the pentagon edge, the magnetic state 
becomes parallel to the y-axis (see the central inset of 
Fig. 1). After switching off the magnetic field, the vor- 
ticity will be clockwise, see the right inset of Fig. [T] To 
describe the vorticity quantitatively, one can introduce 
the chirality similarly to Ref.^: 
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where Xi = cos(27ri/iV) and yi — sm{2Tri/N) if the co- 
ordinate origin is chosen at the center of the pentagon, 
M,^i and Myi are x and y projections of magnetization 
of the i — th particle, respectively. 

To find the area of parameters where the fastest and 
the most reliable reversal occurs, as the characteristic to 
be studied let us choose the first passage time of a certain 
boundary, e.g., when the chirality, starting from the pos- 
itive value, see Fig. [U crosses its -1/2 value level. The 
mean first passage time (the mean reversal time, MKT) 
r, and the standard deviation of the first passage time a 
(SD, jitter) are2 ^: r={t)=j:l,t,/K, {t')=j:tit! / K , 
a—\J (t^) — (t)^, where ti is the first passage time of an 
absorbing boundary and K > 10000 is the number of 
realizations. 

In the calculations^'* let us choose the same param- 
eters as in?-: a=G.l, 7=1.76 • lO^Hz/Oe, Afs=1400 Oe. 
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FIG. 1: The temporal evolutions of magnetization (dashed 
curve) and chirality (solid curve). 



For modeling we take different amplitudes of the mag- 
netic field in the region from i/o=1000 Oe to i/o=2800 
Oe, such that the field 7Jo=1000 Oe slightly exceeds the 
magnetic dipole interaction field. In the following we 
will study different radii of the particles, starting from 
lOnm to 30nm. As an example of a driving with smooth 
fronts we consider the sinusoidal pulse He=eHQ simrt/tp 
with the width tp, where e is the unitary vector of mag- 
netic field direction. If the switching during tp docs not 
happen, the computation is continued for He=0 until 
some long period of time t f , much larger than any other 
relaxation time scale. Our aim is to find the parameter 
range, for which the reversal by a smooth pulse takes 
place quickly and reliably. 
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FIG. 2: The MRT (solid curves) and SD (dashed curves) 
versus pulse width for different angles between the pentagon 
edge and external magnetic field, the temperature r=300 K. 

Let us start our study examining the dependence of 
the MRT and SD versus the angle d between the pen- 
tagon edge and the driving field. Presence of fluctuations 
or small asymmetry of magnetic field will lead to a fast 
formation of the vorticity ([3]), see Fig. [T]for Q — \° and 
T = 300 K, where the initial magnetization was uniform 
as in central inset of Fig. [1] and the vorticity was formed 



in 1 ns. In Fig. [5] the MRT and SD are presented versus 
the pulse duration for different angles between the pen- 
tagon edge and external magnetic field. One can see that 
this dependence is very similar to the single dipole case, 
compare with Ref,^: both MRT and SD demonstrate 
minima as functions of the pulse width and there is a 
strong dependence on the angle 0. The decrease of tem- 
poral characteristics at large tp is due to the fact that 
with decrease of the pulse width the potential barrier 
disappears faster. With further shortening of the pulse, 
the magnetization does not have enough time for the 
complete reversal during ip, so the MRT increases. In 
particular, the reversal time can not be smaller than the 
formation of vorticity, which is of order of 1 ns, see Fig. 
[TJ This, actually, means that for rather short pulses 
the transition occurs due to effect of fluctuations (the 
so-called noise- induced switching). The same applies to 
SD, since for smaller switching time the noise has less 
possibility to move the system far from the deterministic 
trajectory, existing in the absence of noise. The existing 
of the optimal driving angle is again similar to a single 
dipole case and is explained by existence of an unsta- 
ble state of the system, leaving which seriously delays 
the process of switching. One can see that changing the 
driving angle from 9 = 1° (whose value is close to the 
unstable state 6* = 0°) to 6* = 15° allows to speedup the 
switching process by more that one order of magnitude, 
and to correspondingly reduce the SD. However, in dif- 
ference with the single dipole case, the optimal angle 
here is not = 45°^^^, but roughly 6* = 10°, which 
is 1/8 of the angle between two subsequent dipoles in 
the pentagonal array. 
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FIG. 3: The MRT (solid curves) and SD (dashed curves) 
versus the pulse width for different driving pulse amplitude 
Ho, the temperature T=300 K, 9 = 10°. 

An unusual situation in comparison with a single 
dipole case occurs if one tries to study the dependence 
versus the driving field amplitude. This dependence is 
presented in Fig. [3]for the optimal driving angle 9 = 10°. 
Here, larger driving field delays the switching. For ex- 
ample, the increase of the field amplitude from Hq=1AQQ 
Oe to i?o=2800 Oe increases the minimal MRT in almost 
two times. Nevertheless, the minimal attainable SD is 
nearly the same and is proportional to \/T as it must. 



Finally, let us study the dependence from the radius 
of the particle, see Fig. |4l which is very important to un- 
derstand the potential applications of such vortex states 
as components for MRAM. Note that, due to the nor- 
malizations chosen, the only dependence of the noise in- 
tensity on the particle radius remains. As one can see 
from Fig. 21 formally, even for i? = 15 nm for the given 
parameters the reversal still occurs stably. To check this 
more precisely, let us perform the simulation of non- 
switching probability, which for many important appli- 
cations must be below 10~^^. Such low probabilities can 
not be computed numerically due to enormous calcula- 
tion time. Let us perform simulations for 10'°^ realiza- 
tions, and calculate the nonswitching probabilities p{R) 
for different pulse durations as functions of the particle 
radius. The results of simulations are presented in the 
inset of Fig. |4l For A = 0.9 we choose tp=lQ ns and 
tp=2Q ns, while for A = 1.0 - tp=5 ns, and tp=lQ ns. As 
one can see from the inset of Fig. HI the results p{R) for 
different A seriously differ, while the results for different 
tp, but the same A, actually agree. It is clear from the 
presented digits and linear fit of the curves that even for 
A = 1.0 the probability lO"'^^ is reached for i? « 20 nm. 
While the considered parameters range corresponds to 
the regime with significant impact to noise, and below 
10~^ the nonswitching probability can in fact decrease 
faster, it is recommended to use magnetic particles with 
larger values, and to increase the magnetic radius up 
to the particle radius, A — >■ 1, in order to decrease the 
particle radius to 20 nm. 
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FIG. 4: The MRT (solid curves) and SD (dashed curves) ver- 
sus the pulse width for different values of the particle radius 
and A for the temperature T=300 K, 6* = 10°. Inset: the 
nonswitching probability versus particle radius for different 
A and tp 

To understand the reason of rather large non- 
switching probability, let us investigate the distribution 
function 'w{C) of chirality for different values of the par- 
ticle radius and A for the temperature r=300 K. One can 
see from Fig. [S]that the distributions are rather broad 
and the distribution for i? = 15 and A = 1 is more nar- 
row than even for i? = 18 and A = 0.9, while equal mag- 
netic radius would be for R = 16.7 and A = 0.9, which 
is presumably due to stronger dipole interaction in this 



case. Introducing relative width of the distribution func- 
tion c — \J (C-^) — (C)^ / (C) , one can see from the inset 
of Fig. [5] that the increase of A from 0.9 to 1.0 decreases 
the width of the distribution w(C) in nearly two times. 
Therefore, the rather large non-switching probability is 
due to noise-induced instability of the vorticity even in 
the absence of the switching signal. 
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FIG. 5: The distribution function of chirality C for different 
values of the particle radius and A for the temperature r=300 
K. Inset: plot of c = \/ (C) - {Cf/ (C>. 



In the present paper the effect of noise on the pro- 
cess of high-speed remagnetization of vortex state of a 
pentagonal chain of five circular magnetic nanoparticles 
is studied by means of computer simulation of Landau- 
Lifshits model. The mean switching time and the stan- 
dard deviation of the reversal between the vorticities 
have been computed. It has been demonstrated that 
with the reversal by the pulse with sinusoidal shape, the 
optimal pulse duration exists, which minimizes both the 
mean switching time (MST) and the standard deviation 
(SD). Besides, both MST and SD significantly depend 
on the angle between the reversal magnetic field and the 
pentagon edge, and the optimal angle equals 9 ss 10°. 
Also, it is demonstrated that the optimization of the an- 
gle, duration and the amplitude of the driving field leads 
to significant decrease of both MST and SD. In partic- 
ular, for the considered parameters, the MST can be 
decreased from 60 ns to 2-3ns. Such a chain of magnetic 
nanoparticles can effectively be used as an element of 
magnetoresistive memory, and at the temperature 300K 
the stable operation of the element is observed up to 
rather small size of nanoparticles with the raduis of 20 
nm for the chosen experimental parameters. 
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